Application of asymptotic expansions for maximum likelihood estimators' errors to 
gravitational waves from inspiraling binary systems: the network case. 
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This paper describes the most accurate analytical frequentist assessment to date of the un- 
certainties in the estimation of physical parameters from gravitational waves generated by non 
spinning binary systems and Earth-based networks of laser interferometers. The paper quantifies 
how the accuracy in estimating the intrinsic parameters mostly depends on the network signal 
to noise ratio (SNR), but the resolution in the direction of arrival also strongly depends on the 
network geometry. We compare results for 6 different existing and possible global networks and 
two different choices of the parameter space. We show how the fraction of the sky where the one 
sigma angular resolution is below 2 square degrees increases about 3 times when transitioning from 
the Hanford (USA), Livingston (USA) and Cascina (Italy) network to possible 5 sites ones (while 
keeping the network SNR fixed). The technique adopted here is an asymptotic expansion of the 
uncertainties in inverse powers of the SNR where the first order is the inverse Fisher information 
matrix. We show that the commonly employed approach of using a simplified parameter spaces and 
only the Fisher information matrix can largely underestimate the uncertainties (the combined effect 
would lead to a factor 7 for the one sigma sky uncertainty in square degrees at a network SNR of 15). 



I. INTRODUCTION 

With the next generation of advanced ground-based 
gravitational wave detectors under construction, we 
analyse the accuracy in the estimation of physical 
parameters from gravitational waves generated by non 
spinning binary systems. In particular we examine the 
benefits of proposed extensions of the existing LIGO 
Virgo detectors network of Hanford (USA) , Livingston 
(USA) (^38j) and Cascina (Italy) (iilj), to include sites 
in Australia (gD), Japan (gl]) and India ([13]). The 
paper covers the role of the network in the ability (a) to 
estimate intrinsic physical parameters and (b) localize 
sources. This second capability will be critical to aid 
searches for electromagnetic counterparts to detected 
gravitational-wave signals. The analysis is carried out 
with a novel analytical technique developed in [1] and 
extends existing studies (see for example [T], [5], [33] 
[in] J [IS]) [H] and the references therein) as discussed 
below. The two previous applications of this technique 
by the authors involved a single optimally oriented 
interferometer while estimating physical parameters (a) 
from the inspiral phase of the waveform (T) and (b) 
from the phenomenological inspiral-merger-ringdown 
waveforms generated by black hole mergers ([2])- In 
this work we include: (a) the direction of arrival (b) 
the orientation of the orbital plane with respect to the 
observer (c) the use of a network of interferometers 
and (d) parameter spaces of different sizes. This 



aspect needs to be addressed carefully since Fisher 
information matrix can become singular (this is also 
discussed in see [25]). 

Recent frequentist analytical results addressing 
the direction reconstruction, see for example [Ml EH] 
and the references therein, limited themselves to quan- 
tify the accuracy in estimating the direction of arrival 
by using the inverse of the Fisher Information Matrix 
(IFIM) and a reduced parameter space (the angles 
were estimated from the times of arrival assuming the 
intrinsic parameters known). Other articles used the 
full parameter space, but relied on the IFIM to calcu- 
late the errors [SS] [36] . This can be problematic as we 
discussed in [T], [2] and here. In fact, the IFIM delivers 
accurate error predictions in large signal to noise ratios 
(SNRs), however for moderate to low SNRs, where 
the first detections are likely to be recorded, they can 
grossly underestimate the uncertainties. This effect is 
exacerbated when the parameter space is expanded to 
include intrinsic parameters. The role of the size of the 
parameter was initially discussed for IFIM in T2]. 
In [I],[l] we shown how the SNR necessary to trust 
the inverse of the Fisher information matrix strongly 
increases when the parameter space is enlarged to 
include all the unknown parameters in the problem. 
We also discussed how at low SNRs the errors (defined 
as the square root of the sum of the first two order of 
the expansion of the covariance matrix) decreases like 
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1/SNE? this is the SNR region where the the inverse 
of the Fisher information matrix is smaUer than higher 
order contributions we introduced in [1]. 
In the analysis we observe that increasing the size of 
the network from 3 to 5 IFOs can increase by a factor 
3 the fraction of the area that has angular resolution 
of the order of one or few square degrees. We also 
notice that the effects of the size of the parameter 
space and the second order corrections on the direction 
reconstruction can lead to an increase the uncertainties 
by about a factor 5 for a network SNR of about 15. 
The results of section IV show how the accuracy in the 
direction reconstruction strongly increases when the 
network transits from 3 to 4 interferometers but less 
noticeably with further increases. On the positive side 
the details of the analysis show that for regions in the 
sky where the IFIM predict the smallest uncertainties 
the impact of the second order and of the size of the 
parameter space are reduced. 

The asymptotic expansions used in this paper can 
be used to quantify the variance, the bias or the mean 
square error of the maximum likelihood estimators. 
When we talk about "errors" we refer to the square 
root of the mean square error which is identical to the 
square root of the variance (or the diagonal elements of 
the covariance matrix for multivariate estimates) when 
the bias is zero ([!]). Notice that what we call bias is 
not related to errors in the modelling of the signal or 
detector, but simply the shift of the mean value of the 
estimator with respect to the true value of the param- 
eter. The origins of the bias are the non linearity of 
the estimation process and and sometimes edge effects 
(if someone estimates a non negative parameter there is 
usually a positive bias if the true value is close to zero). 
The results of evaluating the expansions of the bias in 
the mass bins considered in this paper show that the 
bias originating from the non linearity of the estima- 
tion process is not a concern, leaving only only edge ef- 
fect biases, when applicable (see |3]). Finally while the 
analysis of this paper is performed within a frequen- 
tist framework, a parallel effort that did not include 
India and Japan studied the problem with numerical 
Bayesian methods ([M])- A direct comparison between 
the two would require to understand the impact of the 
prior probability distributions and the difference in the 
questions the two approaches actually try to answer. 

In section 2 we review the properties of a received 
signal and the conventions in defining the direction of 
arrival and the polarization angle with respect to dif- 
ferent reference frames. The expressions are obtained 
explicitly because none of the many sets of conventions 



adopted in the literature produced suitable expressions 
to be inserted in our asymptotic expansions. In the text 
we provide a more detailed discussion of these issues. 
In section 3 we extend the formalism developed in [1] 
and 2 for a single, optimally oriented interferometer, 
to a network of interferometers. In section 4 we define 
the waveform and the parameter space. In section 5 
we describe the results. In section 6 we provide some 
conclusions and recommendations. In the appendixes 
we provide some supporting proofs and mathematical 
expressions. 



II. THE SIGNAL ON A EARTH BASED 
NETWORK OF DETECTORS 

A. Frames and conventions 

In order to calculate the form of a gravitational wave 
signal in each detector of a network, it is useful to con- 
sider several coordinates frames, as some aspects of the 
problem may become more clear in a frame than in an- 
other. In particular we shall find useful to work in the 
following frames: 

• Wave Frame (WF). It has coordinates {X,Y,Z). 
The GW travels along the positive Z axis; X and 
Y lie in the directions of the polarization ellipse's 
axes. 

• Earth Frame (EF). It has coordinates 
{^EtVEtZe)- The origin is on the center of 
the Earth, the xe axis lies in the line joining the 
center of the Earth with the intersection between 
the Equator and the Greenwich meridian; ze 
passes through the North Pole. 

• Fth-Detector Frame (IDF). It has coordinates 
{xi,yi, zj). The axes xj and yi lie along the two 
arms, n' and n\ respectively, zj points radially 
out from the Earth surface. 

• Fiducial Frame (FF). For the moment is a generic 
frame with respect to whom the positions of the 
detectors will be specified. Later in the section we 
identify the fide with the Earth frame. Henceforth 
we will refer to it as the "fide" frame. 

We introduce two sets of Euler angles: the first one 
transforms the fide frame into the wave frame, the sec- 
ond transforms the fide frame into the I-th detector 
frame: 
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(a^/3^7') 



FF 
FF 



WF 
IDF 



(2.1) 
(2.2) 



Following Goldstein ([H]), and a big part of the lit- 
erature, we perform the second rotation around the in- 
termediary X axis (ZXZ convention): 



(2.3) 



It must be underlined how the Euler angles <j) and 9 
are not to be identified with the spherical coordinates 
of the gravitational wave source in the fide frame. 
If one calls $ and O the spherical coordinates of the 
source in the fide frame, S and ( the angles that the 
wave's direction of propagation forms with the fide axes, 
the following relations hold; 
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(2.4) 



these relations, in terms of the latitude lat. and longi- 
tude long, of the source in the fide frame, become: 



long 



TT 
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+ lat = (. (2.5) 



One can find similar relations between longitude f2^ 
and latitude T-^ of the I-th detector beam splitter and 
the Euler angles , (5^ . The third Euler angle, can 
be fixed noticing that after the second rotation the in- 
termediate frame will have its x axis along the local 
parallel, pointing East, and the y axis along the local 
meridian, pointing North. The third Euler's rotation 
must align the axes with the arms. We label the arms 
in such a way that the first arm has only the x compo- 
nent, while the second has only the y component. If we 
call A-'^ the angle between the first arm and the local 
North direction, then: 



vr 
2 



(2.6) 



The local Earth radius at the position of the I-th 
detector's beam splitter is indicated with rj. 

In the following sections, necessary different usages 
of indexes are introduced as below: 



when the range is 1..3 (Greek letter if the range is 
0..3). The Einstein convention for these indexes 
is always used, unless the contrary is explicitly 
stated. 

2. When a quantity refers to the I-th interferometer, 
we add / on the upper right side of the symbol. 
E.g. 

3. When we calculate explicitly the components of a 
vector or tensor in a particular frame, the frame 
where the expression holds is indicated with an 
index in the upper left of the symbol. E.g. ^^^d'^ 
will be used to give the components of the I-th 
detector tensor in the I-th detector frame. 



B. Interferometer output 



Gravitational waves are ripples in the space-time, 
that can be represented as fluctuations /i^i^ around 
the background flat metric. In the transverse-traceless 
(TT) frame ([S], |3]) the only non vanishing components 
of hp^i,, if the wave is traveling in the positive Z direc- 



tion, are hxx = 
where /i+ and h 



lYY 



Ha. and h 



XY 



l-YX 



are the two independent polarizations. 
While dealing with gravitational waves, it is useful to 
introduce a 3 x 3 wave tensor, defined as one half of the 
spatial part of the metric perturbation: 



i,i = 1..3. 



In an arbitrary frame, the wave tensor can be ex- 
pressed in terms of its circular polarization components 

my- 



{h+{t) + ih: 



,{t))e'^ + {h+{t)-iK 



L 

where r are the right and left circular polarization 



tensors unity, 
to the other: 



They are orthonormal one with respect 



e^e 



= 1, 



^R.^L ij — 



= 1 (2. 



Since the components of the wave tensor are real: 



1. Tensor indexes. Lower case Latin letter, usually 
(i, J, fc, ..), will indicate the components of tensors 



^ Some authors, notably [IS) and |19| . put an explicit imaginary 
unit in the cross polarization, then writing hxY = hyx = ihx 
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(2.9) 



'''/ = ilDF — tpF 



n- (ri - rp) 



(2.17) 



In any reference frame, the polarization tensors can 
be expressed as: 



e'l = rn^rn? (2.10) 

where are the components of the null vector: 

1 
72 



(2.11) 



with ex and unit vectors along the WF axes. Us- 
ing (2.9| and (2.10), the wave tensor (2.7| can be writ- 



ten in an alternative way, that will be useful later: 

w'^t) = hj^(t)W(mimj\ + hy,(t)%mim^\ (2.12) 
In the wave the null vector m has the simple form: 



1 
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(l,z,0) 



(2.13) 



In the previous equation, v\ and are the positions 
of the I-th detector's beam splitter and the fide frame 
origin in an arbitrary frame, while n is a unit vector 
in the direction of the wave propagation, i.e. on the Z 
axis. The time shift is a scalar, and can be calculated 
in any frame. However, the calculations are simpler in 
the fide frame, where rp = and 



PF " ^ 

n(^,() = (sin(^cosS, sinCsin^, cosC) (2-18) 

^^ri(r/,r2^T^) = ri (cos cos sin cos sin T^) 

(the relations between these angles and Euler's were 
given before, eqs. (2.4| and (2.6)). 



From eq. (2.16) it follows that each detector tensor, 
has a very simple expression in its own frame: 



and the wave tensor has the explicit expression: 



If we introduce the null vector: 




(2.19) 



WF„, 







(2.14) 



The form of the wave tensor in another frame is easily 
calculated expressing the components of the null vector 
m in that frame and using eq. (2.12). The GW signal 



at the I-th detector can be written as: 



(2.15) 



where is the I-th detector tensor f|15p. For de- 
tectors with orthogonal arms, like those we consider in 
this work, d has the form: 



(2.16) 



n-*-! and n'2 are the unit vector along the I-th inter- 
ferometer arms. 



In eq. (2.15), the wave tensor is evaluated at the 
retarded time t — tj that takes into account the time 
the GW needs to travel from the I-th detector site to 
the fide frame: 



P =1/2(^-1 



^ e 



yi) 



(2.20) 



where and Gyj are the unit vectors along the IDF 
axes, then the detector tensor can be written in any 
frame as: 



4 =25i[pfp}] 



(2.21) 



just expressing the components of in the frame of 
interest. 



The signal (2.15) can be calculated explicitly in any 
frame. Being an invariant, the final result will be the 
same, but the calculations can be somehow simpler in 
a particular frame. 

Let us consider the wave frame, for example, in which 



the wave tensor has the simple form (2.14). We need to 
write the signal at the I-th detector in the WF, which 
can be done calculating the components of the null vec- 
tor in the WF: 



and using eq. (2.21). The signal acquires the form: 
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Pii Pii 



WF^ WF„ 
P22 P22 



(2.22) 



We do not derive here the exphcit value, as it is rather 
cumbersome. 



Given the general formula, (2.221 we can recover the 



single detector result in the following way: suppose the 
I-th detector is the only one present, so that we can 
identify its frame with the fide frame. This implies that 
the angles and the time delay tj are zero. 



I 

and that the matrix i?^^(a^, /3^, 7^) is the identity ma- 
trix. With some simple algebra is possible to show that 
in this case the signal (2.22 1 can be written as follows: 



s{t) ^ h+{t)F+ + h^{t)F^ 



(2.23) 



where the antenna pattern have the following explicit 
expressions in terms of the Euler angles: 



F+ = - (cos^ 6 + 1) cos 2(/) cos 2ip — cos 6 sin 2(j) sin 2-0 
Fx = — - (cos'^ 6 + 1) cos 2(j) sin 20 — cos 6 sin 20 cos 2-0 



(2.24) 
(2.25) 



The expressions obtained are not immediately com- 
parable with other present in the literature, due to the 
many different conventions used in the years. Let us 
just compare our result with some notable articles. Our 
antenna patterns become the same of [15) if we express 
them in terms of S and C (eq. 2.4), because [15j writes 



them as functions of the angles between the propaga- 
tion vector and the fide axes. More attention is required 
in order to compare our result with those in |18| and 
[19]: in these articles the x axis is not along one of the 
detector's arms, but instead it bisects the arms. This 
introduce a 7r/4 shift in the definition of the first Euler's 
angle. But there is more: in the cited paper the first 
arm is in the positive xy quadrant, while the second arm 
has a negative y component. This means that a relabel 
ni o 7i2 of the arms is also required. This operation 
changes the sign of the detector tensor, and then of the 
antenna patterns. Multiplying eqs. |2.24| and |2.25| by 
— 1, and adding a 7r/4 to one proofs that our antenna 
patterns are equal to those in the references (an extra 
imaginary unit must be added in the cross polarization, 
see footnote [T] at the beginning of this paragraph) 



C. Expansion on a Symmetric- Trace-Free base 

In [18], [19] and [6], the problem of the signal form 
at the detector is treated in an elegant way, recogniz- 
ing that both the detector and the wave tensors are 



r 



symmetric and trace-free (STF) tensors, and can be de- 
veloped in a base of STF tensors. 

We have seen that the wave tensor can be expressed 
in any frame in terms of the real and imaginary part of 
the matrix m m, (the matrix whose (i, j) element is 
minnj). In the wave frame m has the components given 
in eq. (2.13), while in a different frame, say the FF, 



which can be obtained from the WF with rigid rotation, 
its components will depend on the Euler angles that 
describe the rotation. More precisely, it can be shown 
that the real and the imaginary parts of the matrix 
m m in the FF are: 



FF 



[T2„(0,^,V)+T_2„(0,e,0)]:V, 



2J 



(2.26) 



FF 



[T2„ (0,0,7/') -T_2„(0,^, 

(2.27) 

Where the functions T^n are the second-order 
Gel' f and functions, and 3^"^ are a base for rank 2 STF 
tensors. Both these quantities, together with some use- 
ful mathematical background, will be discussed in ap- 
pendix [X] 

As told before, both the wave tensor and the I-th 
detector tensor are STF, and can thus be written in 



terms of the Gcl'fand functions, using (2.12), (2.211 



(2.26) and (2.27) with the appropriate rotation angles. 
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Let us calculate, for example, the signal in the wave 



frame, were the wave tensor is (2.14) and the I-th de- 
tector tensor has the form: 



WF,/ 



/ Ott / 27]- 

25R[p|p}] = 2^— [T2„(WF ^ IDF) + r_2„(WF ^ IDF)]X- ^ 2^Y^XnXJ 



(2.28) 



where we have indicated with the syntax i^(WF — ?> 
IDF the fact that the function F depends on the Euler 
angles that rotate WF into IDF, and where we have 
introduced the combination: 



Xn = T2„(WF ^ IDF) + r-2„(WF IDF). (2.29) 



Using the symmetry (A28) of the Gel'fand functions. 



it is simple to verify that 

* 

X—n — Xn 



(2.30) 



This allows us to write the detector tensor in the WF 
in the simpler form: 



WF, 



dij — 



3?[X2] ~5[X2] 




(2.31) 



Multiplying ( |2.14[ ) by ( |2.31[ ) the signal at the I-th 
detector takes the form: 



/(t) = /l+3?[X2] - /lx3[X2] = 5R[/'^X2] 



(2.32) 



with h \ 



This expression looks beautifully compact, but it is 
not very useful for real calculations. Let us write down 
the explicit value of X2 in terms of Gel'fand functions: 



X2 = r22(WF ^ IDF) + T_22(WF ^ IDF) = 
= r2,(FF ^ IDF) T,2(WF ^ FF) + 
+ r_2s(FF ^ IDF) T,2(WF FF) 
= r,2(WF ^ FF) [T2,(a^ 7^) -fT_2s 7')] 



where we have used the addition formula (A31) for 



Gel'fand functions, while passing from first to second 
line, and the relation (A29) in the last line, to write the 



Gel'fand functions associated with the inverse rotation 
FF ^ WF. 

The real an imaginary part of X2 can be then written: 



J-| = 5R[x2] - 2 (7^2.(a^/3^7') +T-2s(a^/3^7')) {Tlsi^,9,^j) +Tl,MO,i^)) 
Fi = -3[X2] = ^ (T2,(a^/3^7^) +T_2.(a^/3^7')) {T;^{ci^,e,^) ~ Tl^M^.i^)) 



(2.34) 
(2.35) 



where we have introduced the generalized antenna 
pattern and Fy, , using which the signal at the I-th 
detector looks formally as in the single detector case: 

s^{t) = h+{t-Ti)F+ + hy{t-Ti)Fy (2.36) 



when a single detector is present in the Network, this is 
done in appendix [B] The expression we have obtained 
here is compatible with [6] once the differences in the 
definitions of the detector frames are taken into account, 
which results in a 7r/4 shift in the definition of 7^. 



We won't develop these expressions, because of their 
size. Instead we can check that the generalized antenna 
patterns have the single IFO limit values ( 2.24 and 2.25 ) 
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III. FISHER MATRIX, CRLB AND HIGHER 
ORDERS 

We write the output of the I-th detector as the sum 
of the signal {t) and the noise n^{t): 



the sky (see section IV A). However, here and in what 
follows, to unburden the expressions, we do not write 
explicitly this dependence, as well as the dependence on 
the time delay t/, by using the notation: 



(3.1) 



The signal will generally depend on a vector of 
unknown parameters i9 which we want to estimate, that 
can be the physical parameters of the source (e.g. total 
mass) as well as extrinsic parameters as the position on 



s'{t) = s'{t-Tj;^). 



(3.2) 



If we assume that the noise in each detector is sta- 
tionary and Gaussian with zero mean, the probability 
distribution for the data realization in the Tth interfer- 
ometer is given by: 



p{x') cx exp <^ - / dtdt'{x\t) - s\t))VL\t - t'){x'{t') - s'{t')) 



(3.3) 



r 



where Vl^ {t — ti) is the inverse of the noise correlation 
matrix for the I-th detector. The noise from laser in- 
terferometers is usually a combination of smaller Gaus- 
sian fluctuations and larger, rarer non Gaussian outliers 
("glitches" in the data). The use of coincident require- 
ments between different sites and a whole set of data 
quality and vetoes procedures make it reasonable the as- 
sumption that glitches will be recognized an removed or 
vetoed from the data, leaving only Gaussian distortions 
to the GW signal. We also assume that the noise in 
different IFOs is statistically independent. In this case 
the join t probability distribution is simply the product 
of ([33| for / = l..iV: 



N 



P{^) = \[p{x') 



(3.4) 



1=1 



and the log-likelihood turns out to be the sum of the 
log-likelihood of the detectors: 



where E[-] is the expectation with the joint pdf (3.4 1 
f ^[-J is the expectation with the Tth marginal proba- 
bility and where we have defined: 



dd' ■ ■ ■ 3^3 



(3.7) 



The same notation will be used for the signal deriva- 
tives with respect to the components of the vector of 
unknown parameters i?, e.g. Si = 

The single detector Fisher matrix Tj^ can be ex- 
pressed in the Fourier domain and it is possible to prove 
(see for example [T]) that the final expression is: 



(3,8) 



N 



N 



\n[P{^)]^£^Y.^n[p{x')]=J2^' 



1=1 



1=1 



The additivity will hold for functions that are built 
from (3.51 through linear operations. This is true in 



particular for the Fisher information matrix defined as: 



N 



1=1 



N 

Erf. 

1=1 



where we have introduced the I-th detector one-sided 

(3.5) 

noise spectral density {f) ('|16|. |15j). The lower cut- 
off frequency, flow is a detector dependent quantity that 
we fix when the noise models for the detectors are in- 
troduced. The upper limit of integration f^p, which is 
waveform dependent, indicates up to which frequency 
we are confident that the waveform we use is correct. 
These limits will be set later in the section. 

(3.6) Noise-weighted integral like this one, between deriva- 

tives of the signal will, are denote with wedge brackets: 
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{sa-t\sj...p) = {a---i\j---p) = { — 



nfcut 



Sif) 



d6^---de'^d9i ■■■ddP' 
(3.9) 



The optimal SNR, , of the signal at the I-th detec- 
tor can also be written in a similar way (|15j) : 



/\ 2 



45R /"'d/^^^||^ = (.(/)|.(/))^ 

. . (^-^"^ 
It is intuitive that adding interferometers to our net- 
work, while keeping everything else fixed, will increase 
the amount of information we have on a signal that is 
present. We can take this into account introducing the 
network SNR, p defined as: 



N 



(3.11) 



1=1 



which is the SNR we quote in the results section. 
Once the Fisher matrix is calculated, its inverse can be 
used to estimate errors and covariances. 

The Cramer Rao lower bound (CRLB) gives a bound 
for the variance of the i-th parameter f|17|): 



(Ai^*)2 > [r-i]^^ = r" 



(3.12) 



while the covariance between the i-th and the j-th 
parameter is: 



(3.13) 



In a high SNR regime, the errors and covariances are 
close to these bounds, and one writes the probability 
distribution for the errors as a multivariate Gaussian: 



p(Ai9) cx 



(3.14) 



with i,j = 1..M, M being the dimension of i9. 

On the other hand, the CRLB is known to fail for 
small SNRs, that is, for weak signals. In those situ- 
ations numerical simulations, like MonteCarlo simula- 
tions (MC) are more faithful, but at the price of an 
high computational cost. In [1] we proposed an analyt- 
ical method to improve the errors estimation for small 
SNRs, and we used it for a 3.5PN signal. In [2] we ap- 
plied the same method to an Inspiral-Merger-Ringdown 
(IMR) signal ([9-11]), also showing that the bias arising 
from the nonlinearity of the estimator plays a funda- 
mental role for high mass systems. 

In ^ and [2] we have shown how the bias and 
(co)variance of the estimators can be written as power 
series in one over the SNR p: 



4W 

B^[l] 



p2 



6^.[l]+&^.[2]-|-...(3.16) 



where after the first equal sign we have shown explic- 
itly how the different terms depend on the SNR. We 
have also shown that the first order in the variance se- 
ries is the usual CRLB, while the second orders contain 
higher derivatives of the signal. Formally, we can still 
use the expressions we gave there, for example for the 
second order diagonal elements of the covariance ma- 
trix: 



_|_ pjmpjnpP^pqt 
+ Vnq,z 



(3.17) 



but now both the network Fisher matrix inverse F*^ likelihood, 
and the v must be calculate using the joint network 

The second corrections contain the factor: 
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where now i is the jomt log hkehhood, eq. (3.5). It 



is easy to check that it is not possible to just add up 
the single-IFO contributions, as "cross terms" between 
different interferometers might appear. 
Let's consider for example: 



The terms with I ~ J will give back the sum of single- 
IFO f^f, but the terms with I ^ J also give a con- 
tribution, for example for the (1 — 2) term (note the 
apex at the end of the square brackets to distinguish 
the contributions of IPO 1 and IFO 2): 



N N 



Vab,cd 



= E [LbU =Y.J2^ KbC] ■ (3-18) 



/ J 



Sab{t, - t') [ - S{t', d) ) - Salt, §)n{t - t')sb{t', d) 

1(2) 



(1) 



E Kb^ld] = j dMi'dfcdfc'dxP(x) 

X s^a{s, ??)f^(fc - k') ^x(fc') - s{k', t?) j - Sc{k, d)n{k - k')sa{k', d) 

= f dtdt'dx^pix^)si^\t,^)n'^^\t-t')si^\t',^) f dkdk'&x''p{x'')sf\kj)n^^\k~k')sf{k'j) 



— ^ ab^ cd 



(3.19) 



It is then clear that (3.181 can be written as 



N 



N 



N 



Vab,cd = E <b,cd + E ^ab^cd (3-20) 
7=1 iTtJ 



Similar calculations give the following results: 



N 



Va,b = -t^a 



= Tab ~ y^r^fa 

1=1 

Vab,c = ^ , c 
1=1 
N 

E ^'abc , d 



N 



Vabc , d 



Vabc 



1=1 
N 



abc 



(3.21) 
(3.22) 
(3.23) 
(3.24) 



1=1 



Vab , cd 



Vabcd 



Vab .cd 



^abc , de 



Vabcd , 



Vabc , d , e 



Vab , cd , 



Vabcde 



E ^"-b . ^ E '>^' 
1=1 I^J 

N 

E-i 



J 

c d 



abed 



1=1 

N 



E'^af'.c.rf E ^ab^cd 



1=1 



N 



abc , de 



abed , 



J 

abc 



E-i 

1=1 

N 

E-i 

1=1 

N 

E ^^f'c ,<>.,e + E ^abc^de 



1=1 
N 



1^.1 



(3.25) 
(3.26) 

(3.27) 
(3.28) 
(3.29) 

(3.30) 
(3.31) 



E-i 

1=1 

N 



ab , cd , e 



abcde 



J2Hb,c^id + yUe^P) 

(3.33) 



1=1 
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where the single-detector U| j have the values 
in [2], that we report in appendix [P] 


given 












N 




The V in which cross terms are present can be further 


Va,b 


= -Vab = Tafc ^ r^fc 

T 1 

1 i 


(3.38) 


simplified. Let us consider again Vab,cd- 






N 








Vab,c 


= ^ ^^b , c 
7=1 


(3.39) 


N 

Vab , cd = ^ Uab , cd + ^ b^cd 


(3.34) 


Vabc , d 


Af 


(3.40) 


1=1 I^J 






7=1 

AT 




The single detector has the value: 




Vabc 


- v,/ 

— ^ '^abc 








Vab , cd 


= 7 (Sab , Scd)"'^ + Tabbed 


(3.42) 


''^ab.cd = {^ab , ScdY + ^ab^id- 


(3.35) 


Vabcd 


7=1 

Af 

— ^ ^abcd 


(3.43) 


SO that eq. |3.34|can be expressed as: 




Vab ,c,d 


7=1 
= —Fa bred 

Af 

— ^ ^('^gbc 1 -^de) '^abcTde 
7=1 
AT 

/ ^ abed , e 


(3.44) 


N N 
Vab cd — / ^ (Sab , Scd) ^ + TLfL + Vi 
1=1 1=1 I^J 


cd 

(3.36) 


Vabc , de 
Vabcd , e 


(3.45) 
(3.46) 


It is nearly evident that: 




Vabc , d , e 




(3.47) 






Vab cd c 


— ^^-^ db'^ cd £ cd'^ ab e 
N 

/ J abcde' 
7=1 


(3.48) 


AT 

'Y^^ab^cd + ^ ^ab^cd — ^ ab^ cd , 




Vabcde 


(3.49) 



The fact that we have been able to prove that the f 's 
explicitly contain the network Fisher matrix is impor- 
which implies that only the network Fisher informa- tant, as it allows to obtain simplified expressions for the 
tion appears in the final form: first order variance and bias: 



^ a^„[l]=r- (3.50) 

Vab, cd = 'Y{Sab, ScdY +Vab^cd (3.37) 
7=1 

bs^ [1] = ^r™F''^(i;,,e + 2vc,ab) (3.51) 

The same kind of calculations show that no cross- 
terms are present in any of the v, whose final form is: the second order variance: 
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N 



_j_ pj? pj? pP p"? H~ '^Vnpq'^mzt ~^ '^'^qz .n'^mtp ~^ '^'^qpjz'^nmt H" 

^^mqp^nt , z ~^ '^pqz'^nt , m ~f~ '^'^rnq , z'^pt , n ~l~ '^'^pt , z'^mq . n ~l~ ^mz , t'^nq . p 



(3.52) 



and the second order bias and covariance matrix, that 
we show in appendix [C] 



with Tl_ and defined in eqs. (2.34 1 and (2.351. 
The Fourier transform of the signal considering 
only positive frequencies, is: 



IV. RESULTS 
A. Gravitational waveforms 

In this paper we analyse the errors in the estimation 
of the physical parameters of a G W signals generated by 
the inspiral phase of binary systems whose components 
are not rapidly spinning. We assume that when the 
frequency of signal enters in the bandwidth of the ad- 
vanced detectors, these systems have already lost their 
orbit eccentricity. In this case the plus and cross polar- 
izations are (gj, [15]): 



1 



h{t)+ = A{t) 



h{t)x = A(t) cose sm^{t) 



cos^{t) 



(4.1) 
(4.2) 



where A{t) is a time dependent amplitude, $(*) is the 
phase. Both these quantities can be calculated with a 
higher degree of precision within the frame of the Post- 
Newtonian theory ([5], |2lJ- The angle e is defined 
later in the text. The signal (2.36) can be, in this case, 
written in the form: 



s\t) = / A{t - Ti) cos($(i - Ti) - $^) (4.3) 



where 



1 + cos^ e , 



<I>g = arctan 



2 cos e 
1 + cos2 e)Fi 



{coseFlf (4.4) 



(4.5) 



dtAit~Ti)er^-'''^'+"^^'-^'^~'< = 
-■iTTifn-i-ii I dT^(T)e^^''^'^-''+**'^^ (4.6) 



The integral can be solved using the well known sta- 
tionary phase approximation - SPA - ([12], [21]) which 
consists in expanding the integrand around its station- 
ary point, where the derivative of the phase is zero. 
This also explains why we only kept e~*^*'-*~'^^^~*o^ 
while going from the first to the second line: the second 
term would results in an second integral whose inte- 
grand would have a phase: 



i~27r ft - $(^ - r) + $o) 



(4.7) 



But the derivative of this expression with respect to 
the time is never zero, because the derivative of the 
orbital phase ^{t) is always positive. This implies that 
the phase of this additional addend would be oscillating 
in the whole range of integration, making the integral 
small. The final result is then 



TTiD\ 24 

where the phase is given at the 3.5 PN order by: 



^(/) = 27r/t, + 0,-^ + ^,^^afet;'= (4.9) 



4 128t]v^ 



fc=0 
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and V — {'KMf)3. The coefficients a^, that depend 
on the total and symmetrized mass, can be found in |12| 
and [ig. 

Let us summarize the unknown parameters on which 
this waveform depends: 

• tc and (j>c are a reference time (usually the detec- 
tion time, or the coalescence time) and the phase 
the wave had at that time 

• M = rni + TO2 is the total mass, ri — "11^2 jg 

yTTi 1 -p 77T2 ) 

the symmetric mass ratio. The chirp mass A4 — 
rj^ M is often used instead of the total mass. 

• D is the luminosity distance of the system. 

• t/i is the polarization angle 

• e is the angle formed between the line of sight 
and the system orbital angular momentum (often 
referred as iota in the literature). A system with 
e = is called face-on; one with e — ti/2 edge-on. 
Note that, because of eq. (4.2 ), an edge-on system 



only has the plus polarization. 

Q and $ are the spherical coordinates of the 
sources in the fide frame. It is worth recalling that 
T/ , fj, and $0 are functions of those angles angles 
through (pT7l),(l234| , (plsl) , (lill) and (|45). 



This waveform is usually assumed to be accurate till 
the innermost stable circular orbit (ISCO) frequency. 
This will also be the upper limit of integration in the 
noise-weighted integrals like p. 



ISCO 



fup — 



(4.10) 



It is know, and easy to verify, that the distance and 
e are strongly correlated. From a numerical point of 
view this makes the inversion of the Fisher informa- 



being close to zero. This is a known problem (|29j) 
with explicit or emerging degeneracies. Numerical is- 
sues need to be handled carefully since the second or- 



der of the covariance expansion (3.13) contain the IFIM 
multiplied by itself up to four times. To prevent nu- 
merical instabilities, in this paper we consider two pa- 
rameters spaces: one where only the position of the 
source is unknown (mostly to compare our results with 
existing literature that adopt this simplification), and 
1? = {ta,logAi,ri,lat,long). We checked that in order 
to extend the parameter space further (for example up 
to 9 parameters) the numerical precision would need to 
exceed one hundred digits. 

B. Detectors locations and noise models 

In tablelljwe report the position of the detectors beam 
splitter [23[ with respect to the Earth Frame (eq. (2.2) 
can be used to convert those in the Euler angles needed 
to calculate the generalized antenna patterns). More 
details on the existing and planned IFOs can be found 
in ([3ai3H!,[39H41], [42], [43], [44]). 



Detector 


Label 


Longitude 


Latitude 


A 


LIGO Livingston 


L 


90°46'27.3" W 


30°33'46.4' 


N 


108°0'0" 


LIGO Hanford 


H 


119°24'27.6" W 


46°27'18.5' 


N 


36.8°0'0" 


Virgo, Italy 


V 


10°30'16"E 


43°37'53" 


N 


341°30'0" 


LCGT, Japan 


J 


137°10'48" E 


36°15'00" 


N 


295°0'0" 


AIGO, Australia 


A 


115°42'51"E 


31°21'29" 


S 


270°0'0" 


INDIGO, India 


I 


74°02'59"E 


19°05'47" 


S 


45°0'0" 



Table I. The positions and orientations of the detector in the 
Earth Frame. The angle A is defined as the angle between 
the local North and the first arm of the interferometer mea- 



tion matrix (3.6) numerically unstable, the determinant 



sured counter clockwise. Eq. (2.61 can be used to convert 
those angles in the Euler angles that transform the Detector 
Frame in the Earth Frame. 



For the advanced Virgo detector we have used the 
noise power spectral density plotted in Fig. [ijand given 
in 12j: 



Shif) = So 



2.67 10-^ 0.59 e('"^)' [-3.2-i.o8inx-o.i3(i„.)^]^-4.i ^ O^gg ^-0.73 (in.)=^. 



5.34 



Shif)^^, f<flo 



•) f ^ flow 

(4.11) 



Where x=j-^,fQ^ 720Hz, and Sq = lO^^^Hz^^ 
or the advanced LIGO detector, as well as for the 



projected new IFOs, we use the following noise power 
spectral density, also plotted in Fig. [l] 
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ShW = So 



-4.14 



5x-^ + ni ' 



1 + a;2/2 



1 f ^ flow 

(4.12) 



Where x = 



X EE /o = 215Hz, and 5o = 10-49Hz~\ 
For both the noise spectral densities, we have chosen 
the lower frequency cutoff to be flow = 20Hz. 

Obviously the noise for the new detectors need not to 
be exactly like this one, but using the AdvLIGO con- 
figuration is a good approximation for the sensitivities 
of those new detectors, and the differences will be very 
small compared to the addition of a new IFO to the 
network. 
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Figure 1. (Color online) Square root of the noise PSD for 
an Advanced Ligo detector (continuous line) an Advanced 
Virgo detector (dotted line) 



C. Representative examples and figures of merit 

In the frcquentist formalism the errors depend on the 
value of the parameters. Therefore it is informative to 
show which kind of errors, calculated using the first two 
terms of eq. (3.15), we can expect for different repre- 



sentative signals and different locations in the sky. It is 
also important to keep in mind that some situations are 
intrinsically pathological and less interesting: for exam- 
ple if the latitude is or tt the longitude is undefined 
so it is not meaningful to compute the errors on it (the 
FIM would actually be singular). In thi s paper we do 
not quote results for the bias, eq. (3.16), as we shown 



in [2] that those are important only for systems much 
more massive than the ones we consider in this work. 

The plots in this paper refer to the following binary 
systems: 

• BNS - Binary neutron star - toi — m2 — \AMq 

• BHNS - Black hole - neutron star - mi = 
lOM0,m2 = \AMq 

• BBH - Binary black hole - m\ = m2 = lOAf© 
and network configurations: 

• 3 Detectors: HLV 

• 4 Detectors: HLVJ, HLVA and HLVI 

• 5 Detectors: HLVJA and HLVJI 

When we calculate the errors, unlike some existing 
studies (e.g. [28]), we do not fix the distance of the 
sources to be the same for all the networks. The rea- 
son is that such approach does not disentangle the role 
of SNR gains and network geometry (see some related 
discussion in l54]). Being more interested in the role of 
the network geometry, for each network and mass bin 
we set the distance of the source so that the average 
network SNR over the sky is about 15. In particular: 

• 3 Detectors: BNS: 200Mpc, BHNS: 450Mpc, 
BBH:1100Mpc 

• 4 Detectors: BNS: 250Mpc, BHNS: 550Mpc, 
BBH:1300Mpc 

• 5 Detectors: BNS: 280Mpc, BHNS: 6250Mpc, 
BBH:1500Mpc 

Note however how the the joint effect of additional 
IFOs and higher SNR can be quantified from the results 
we quote in this paper. As we provide both a[l] and 
cr[2], as well as the network SNR, p* , for which they are 
calculated, one can build for any SNR, p, the quantity 
(which is easily derived from eq. 3.15): 



(4.13) 



which gives the errors for the desired SNR. 
In order to study the sky dependence of the angular 
resolution we adopt a 289 points sky grid where the 



14 



(i,j)-tli point has latitude and longitude in the Earth 
frame. 



latj 



long^. 



2 



25 



] 



16 

^ 27r - 25 

TT + d H — — I 

16 



withi,j = 0,1,--- ,16; (5 = 0.02. 



(4.14) 



In each point the Fisher matrix and the second or- 
der corrections are evaluated. The small offset 5 in the 
range of the latitude is added to avoid the poles. The 
first and last values of the latitude are then f — (5 and 

In our simulation we set to representative values of 
the orientation angle e = 7r/6 and ip = Q. Even if the 
results are almost independent of -0 the choice of the 
values of e deserve some discussion. It can be shown, 
from the derivatives of |4.4| and |4.5| with respect to e, 
that for e = derivative of the signal with respect to e is 



zero and therefore the Fisher matrix is singular. Similar 
inversion problems are present for small but non zero 
values of e. However, apart from the neighbourhood 
of zero, the results would not be qualitatively different 
choosing different values of e. We explicitly varied e, 
keeping the other parameters fixed, and observed only 
changes in SNR due to the fact that the cross polariza- 
tion become less and less important (4.2 1. By rescaling 



the distance to compensate for this loss, the pattern on 
the errors in the sky-map are roughly unchanged. This 
is also the case with the larger networks we considered, 
where the polarization angle can be partially resolved, 
and the pattern on the errors in the sky-map is pretty 
much unchanged. 

In the following we describe how the accuracies in 
the estimation of the direction of arrival are be pre- 
sented. The covariance matrix and its second order give 
us the errors in latitude, longitude, and their correla- 
tion. Assuming that the likelihood is roughly ellipsoidal 
the probability distribution of the angular position is a 
bivariatc Gaussian that can be written as: 



p{Mat, Along) cx e-3(^'''*'r,„t,,„t+A/o«g=^r,„„.,,,„„,+2A;atAionsr,„t,,„„g) 



Where F is the inverse of the covariance matrix (or 
FIM if only the first order is being taken into account). 
The contours of constant probability are ellipses in the 
(long,lat) plane. Because of the correlation term, the 
axis of the ellipse are not parallel to the coordinate axes 
but form an angle a. 

One can calculate the principal directions of the el- 
lipses diagonalizing the matrix: 



^ long ^long ^lat,long 
, ^latjlong ^lat.lat , 



(4.15) 



and finding eigenvalues Ai and A2. It is evident that 
the widths of the ellipse along the two eigenvectors are 
then simply 5i = :^j=) with i = 1,2. The solid angle 
corresponding to one of such ellipses, centered around 
a point of latitude 7 can be written as: 



An = 27r| COS7I 



1 



•\A1A2 



(4.16) 



Considering the scalar product of the first eigenvector 
with the unit vector of the long axis one can find the 



angle a. These ellipses are then plotted using an Aitoff 
projection of the sky. 



D. Bidimensional Parameter Space 

In this section we compute the uncertainties for the 
location of a binary system, while the other parameters 
are assumed known. The approximation of reducing 
the parameter space, used for example in [27], [28], [26] . 
unfortunately tends to produce optimistic predictions 
as discussed in details in the next section. The perfor- 
mance of the different network configurations consid- 
ered are summarized in tables [TTl |III| and |IV[ 

While the entries that use just the IFIM in table II 
and for the HLVA network in table III are consistent 
with the results in [2S] , [35] , columns II IV and VI show 
that in realistic SNR regimes the second order needs 
to be taken into account. For example, if we include 
the second order for an average network SNR of 15, 
the solid angle uncertainty increases on average 45%, 
or more in detail between 37% and 137% with 3 IFOs, 
between 10% and 164% with 4 IFOs and between 10% 
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and 153% with 5 IFOs . 

These tables indicate the benefits of enlarging the 
IFO network however we defer interpretation and dis- 
cussions to the next section where intrinsic parameters 
are also estimated at the same time as the direction of 
arrival. 







BNS 


BHNS 


BBH 






a[l] 


a[2] 




a[2] 


a[l] 


a[2] 




lat [mrad] 


7.69 


6.18 


11.3 


6.01 


15.3 


6.85 




long [mrad] 


15.1 


8.13 


20.3 


87.92 


25.6 


8.86 




An[deg^] 


0.84 


2.23* 


1.55 


2.99* 


2.80 


4.45* 




Net. SNR 


16.55 


15.37 


14.26 



Table II. First order, cr[l], and second order, (t[2], errors 
averaged over the sky using a three detectors network. Only 
the sky position is considered unknown 
* Note: for the solid angle error, the column a[2] takes into 
account both the first and second order errors. The reason 



is that for the solid angle error, eq. (4.161, is not possible 



to separate the effects of first and second order. 







BNS 


BHNS 


BBH 






a[l] 


a[2] 


a[l] 


a[2] 


a[l] 


a[2] 




lat [mrad] 


7.96 


6.09 


10.9 


5.02 


14.3 


5.53 


1 


long [mrad] 


14.2 


7.31 


18.1 


6.38 


21.9 


6.82 




An[deg'^] 


0.89 


1.65* 


1.47 


2.11* 


2.47 


3.19* 




Net. SNR 


15.7 


14.94 


14.29 




lat [mrad] 


7.24 


4.12 


9.33 


3.30 


12.3 


3.85 


long [mrad] 


14.7 


6.16 


17.6 


5.03 


21.0 


5.65 




An[deg'^] 


0.77 


1.06* 


1.17 


1.36* 


1.99 


2.23* 




Net. SNR 


15.4 


14.94 


14.05 




lat [mrad] 


6.87 


5.42 


9.09 


4.56 


11.9 


5.00 




long [mrad] 


13.5 


7.39 


16.7 


6.28 


20.4 


6.82 




An[deg^] 


0.75 


2.17* 


1.18 


2.30* 


2.02 


3.04* 




Net. SNR 


16.01 


15.23 


14.57 



Table III. Same as table [III but using 4-IFO networks 



E. Five dimensional parameter space 

In this section we analyze the errors when a five di- 
mensional parameter space i? = {ta,logJ^,r], lat, long) 
is used and discuss in details the benefits of enlarging 
the HLV network of interferometers. Quantitative re- 
sults are presented in tables (|V VI and VII ) and in fig- 







BNS 


BHNS 


BBH 






a[l] 


a[2] 


a[l] 


a[2] 


a[l] 


a[2] 




lat [mrad] 


6.44 


3.46 


7.97 


2.66 


10.7 


3.31 


> 


long [mrad] 


12.9 


4.86 


15.1 


4.08 


18.5 


4.86 


K 


An[deg'^] 


0.65 


0.85* 


0.94 


1.05* 


1.67 


1.84* 




Net. SNR 


15.92 


14.97 


14.06 


< 


lat [mrad] 


6.76 


4.54 


8.72 


3.53 


11.7 


4.24 




long [mrad] 


13.1 


6.04 


16.0 


4.96 


19.7 


5.80 




An[deg^] 


0.71 


1.12* 


1.08 


1.36* 


1.92 


2.29* 




Net. SNR 


16.18 


15.18 


14.28 



Table IV. Same as table IH] but using 5-IFO networks 



ures 2 to 10. The discussion that follow is organized in 
two lines: (a) trends in the accuracy of the direction re- 
construction and estimation of intrinsic parameters (b) 
methodology recommendations to obtain reliable esti- 
mates. 

(a) It is known that with a three detectors network 
the source sky localization is generally not very good, 
and that several blind spots exist in the projection on 
the sky of the plane containing the three IFO. This 
plane corresponds to the low-SNR region of fig. [2] 




Figure 2. (Color online) Network SNR with an HLV network 
for a BNS system at 200Mpc of distance. 

In fig. [3] we plot the ellipses corresponding to 95% 
confidence interval for the BNS system; the solid el- 
lipses are obtained using the first order errors, while 
the dashed ellipses take into account the second order. 
We have plotted a filled contour of the network SNR 
in the background as comparison. It is evident how, 
even though most of the biggest ellipses are in low-SNR 
regions, as someone might expect, there is no common 
proportionality between the two, and the contribution 
of the second order to the sky localization can be big 
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for regions with medium SNR (e.g. long 45 lat 45 S in 
fig. [2). We return to this point later in the section. In 
table V] we report the errors, averaged over the whole 
sky, for the three binary system considered here, when 
detected with the HLV network. 




25.0 
22.5 
20.0 
17.5 
15.0 
12.5 
10.0 
5 



Figure 3. (Color online) Skymap of the errors for a BNS 
signal detected in HLV. The Network SNR is plotted in the 
background and its value is shown in the colorbar 







BNS 


BHNS 


BBH 






a[l] 


a[2] 


a[l] 


a [2] 


a[l] 


a [2] 




M [%] 


4.3Ie-3 


2.58e-3 


2.66e-2 


1.68e-2 


0.21 


0.14 




r, [%] 


0.22 


0.14 


0.38 


0.21 


1.61 


1.00 




ta [ms] 


0.34 


0.23 


0.46 


0.27 


0.68 


0.39 




lat [mrad] 


22.5 


17.9 


28.6 


19.7 


36.1 


22.0 




long [mrad] 


51.9 


46.1 


57.0 


35.3 


76.1 


48.7 






4.71 


16.4* 


9.07 


21.5* 


15.14 


27.09* 




Net. SNR 


16.55 


15.37 


14.27 



Table V. First order, cr[l], and second order, (t[2], errors 
averaged over the sky using a three detectors network and 
a five dimensional parameter space 

* Note: for the solid angle error, the column a [2] takes into 
accout both the first and second order errors. The reason is 



that for the solid angle error, eq. (4.161, is not possible to 



separate the effects of first and second order. 

The sky-average accuracy, at an average network 
SNR of about 15, averaged between the three mass 
bins that we consider here, is 3.2 square degrees for 
HLV when only the two dimensional parameter space 
is assumed. However for the same network and the five 
dimensional parameter space this number raises to 21.6 
square degrees. 

Increasing the size of the network to 4 IFOs improve 
the situation substantially with these two numbers be- 
coming 2.1 and 9.0 square degrees (AIGO is the best 



addition with 5.5 square degree for the 5 dimensional 
parameter space). Finally if the network includes 5 in- 
terferometers these uncertainties would decrease to 1.4 
and 3.5 square degrees. 

It is also important to notice that the angular resolu- 
tion is not homogeneous in the sky (especially with the 
HLV network) . and that the role of the second order in 
the expansions is usually milder in the locations where 
the accuracy is of the order of a degree. 

One way to quantify the fraction of the sky that can 
provide sufficient directional accuracy for electro mag- 
netic follow-ups is described in table |VIII[ For example 
we can observe how the fraction of the sky (averaged 
over the 3 binary signals) where the one sigma angular 
resolution is below 2 square degrees is 20% for the HLV 
network and increases to 62% for the best performing 5 
IFOs network (HLV J A). 

Comparing the results of the four site networks, sum- 
marized in table |Vl] with the three IFO network, table 
|V| it is evident that there are not substantial improve- 
ments in the estimation of the intrinsic parameters, 
chirp mass and t], while the errors in the estimation 
of the arrival time, and sky position are dramatically 
lower. 

Looking at the plots we can also observe that the 
low sky-error regions do not overlap with high-SNR re- 
gions. The reason is that the direction uncertainties 
mostly depend on the geometry of the network and the 
estimation of the arrival time. In fig. |4] we plot the 
error ellipses corresponding to a BNS signal detected in 
HLVA, but this time we put the error on the arrival time 
(first plus second order) as a filled contour background, 
instead of the network SNR. The match between large 
errors in time and position is evident. This also implies 
that if the exact position and orientation of new detec- 
tors is chosen to guarantee the highest network SNR, 
on average, it might not need not to assure the best 
performances for sky localization ([H]). This fact is 
confirmed by table |Vl[ from which is clear that HLV J is 
the network that guarantee the highest averaged SNR, 
but not the lowest errors for intrinsic parameters. 

In fig. [5]and[6]we show errors on the reconstruction of 
the same source, but using HLVI and HLVJ networks. 

Those plots seems to suggest that the HLVA network 
allows the best sky localization, which is also the out- 
come of table VI This is reasonable given that 1) the 
AIGO site is the furthest from the HLV IFOs, and 2) 
the plane containing AIGO's arms would be nearly per- 
pendicular to the HLV plane (providing sensitivity in 
most of the blind spots of the HLV network) . Note that 
HLVA is also the network for which the second order 
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Figure 4. (Color online) Skymap for a BNS signal detected 
in HLVA. The error on the arrival time is plotted in the 
background and its value is shown in the colorbar (in ms) 



75° 




-75" 



Figure 5. (Color online) Skymap for a BNS signal detected 
in HLVI. The error on the arrival time is plotted in the 
background and its value is shown in the colorbar (in ms) 



Figure 6. (Color online) Skymap for a BNS signal detected 
in HLVJ. The error on the arrival time is plotted in the 
background and its value is shown in the colorbar (in ms) 




Figure 7. (Color online) Ratio between the second and first 
order errors for the chirp mass. BNS detected in HLVA. 
Network SNR in background and colormap 



contributions for Afi are smallest. On average the ra- 
tfo between An[lst + 2nd] and An[lst] is 3.15 for BNS, 
2.2 for BHNS and 1.79 for BBH, using HLVA. The same 
ratios are respectively 5.17, 2.76 and 1.92 using HLVI 
and 6.76, 2.95 and 2.13 using HLVJ For the intrinsic 
parameters, chirp mass and r], the second order contri- 
butions are larger were the SNR is lower; in fig. [Tj for 
example, we plot the ratio ^^jfj in contours, with the 
network SNR in the background and colormap. For the 
extrinsic parameters, on the other hand, the second or- 
der terms are larger where the error on the first order 
arrival time were large. This is well shown in fig. |8j 
where it is visible that the biggest contributions to the 
second order errors in latitude come from regions where 
the arrival time error is large. 

The same trends are visible in the two five detec- 
tor networks we have considered: HLVJA and HLVJI. 



The results are summarize in table [VII} Here again it is 
visible how the intrinsic parameters do not take any ad- 
vantage from the addition of a new IFO to the network, 
keeping a comparable SNR. 

The angular errors quoted in table |VII| seem to sug- 
gest that HLVJI is the best 5 IFOs network. However 
one has to remember that what is given in table VII is 



an average of Ail over the 189 sky positions we con- 
sidered. The reason why the first + second order error 
Af2[l -I- 2] is slightly larger for HLVJA is that there are 
a few spots in the sky (four, corresponding to positions 
where the time-resolution is poor) with very large sky 
errors, ~ 60 deg^, that increase the sky-average. These 
spots are clearly visible in Fig. [9j Aside from these 
points, HLVJA is better performing. That is seen in ta- 
bic |VIII[ in which the number of signals with sky error 
smaller than 1, 2 or 3 square degrees are given. 

It is worth mentioning that this trends in the sky 
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Figure 8. (Color online) Ratio between the second and first 
order errors for the latitude estimation. BNS detected in 
HLVA. Error in arrival time (ms) in background and col- 
ormap 



a 4-IFO to a 5-IFO network, the gain in sky localization 
is smaller than going from a 3-IFO to a 4-IFO network. 

In this paper we also observe that the increment 
in performance depend on which detectors have been 
added, and in which order. For example, using only the 
CRLB only and going from HLV to HLVJ to HLVJA, 
°'2o comparing tables |V VI and VII it can be seen that the 
015 first order A17[l] for BNS systems varies in the foUow- 
0.10 ing way (in parenthesis, the variation with respect to 
0.05 HLV): 

• HLV: 4.71 deg^ 
HLVJ: 1.98 deg^ (-58%) 

,2 



localization accuracy, while adding IFOs to the network, 
was investigated with Monte Carlo simulations for burst 
signals in [46]. Explicitly it was noted that going from 







BNS 


BHNS 


BBH 






a[l] 


a[2] 


a[l] 


a [2] 


o[l] 


a[2] 




M [%] 


4.56e-3 


2.92e-3 


2.79e-2 


1.86e-2 


0.21 


0.15 


1 


V [%] 


0.22 


0.14 


0.38 


0.22 


1.58 


0.99 


ta [ms] 


0.11 


0.11 


0.17 


0.12 


0.36 


0.24 




lat [mrad] 


9.58 


9.51 


13.5 


9.78 


17.8 


10.1 




long [mrad] 


39.0 


40.0 


46.6 


31.7 


61.3 


41.5 






1.34 


4.22* 


2.32 


5.11* 


4.01 


7.20* 




Net. SNR 


15.7 


14.9 


14.29 




M [%] 


4.63e-3 


3.00e-3 


2.81e-2 


1.88e-2 


0.21 


0.15 




V [%] 


0.23 


0.16 


0.39 


0.23 


1.61 


1.03 




ta [ms] 


0.17 


0.17 


0.27 


0.18 


0.47 


0.31 




lat [mrad] 


14.6 


15.3 


20.4 


14.1 


26.3 


15.0 




long [mrad] 


46.0 


48.3 


51.2 


34.6 


66.9 


47.4 






1.85 


9.57* 


3.09 


8.50* 


5.27 


10.12* 




Not. SNR 


15.4 


14.7 


14.05 




M [%] 


4.39e-3 


2.64e-3 


2.67e-2 


1.68e-2 


0.20 


0.13 




n[%\ 


0.21 


0.13 


0.37 


0.20 


1.50 


0.89 




ta [ms] 


0.17 


0.15 


0.26 


0.17 


0.44 


0.27 


III 


lat [mrad] 


13.6 


16.3 


19.5 


16.1 


25.0 


17.0 




long [mrad] 


37.4 


39.3 


44.9 


31.7 


58.8 


40.7 






1.98 


13.39* 


3.49 


10.29* 


5.94 


12.68* 




Net. SNR 


16.01 


15.23 


14.57 



Table VI. Same as table |V] but using 4-IFO networks 



• HLVJA: 1.02 deg^ (-78%) 

On the other hand if the order is HLV, HLVA, HLVJA 
we find: 

• HLV: 4.71 deg^ 

• HLVA: 1.34 deg^ (-72%) 

• HLVJA: 1.02 deg^ (-78%) 

The situation is even more clear when the second or- 
der is taken into account. For example if LCGT is the 
first IFO to be added to HLV the errors Af][l -f- 2] vary 
as: 

• HLV: 16.4 deg^ 

• HLVJ: 13.4 deg^ (-18%) 

• HLVJA: 2.78 deg^ (-83%) 

but if the order is HLV, HLVA, HLVJA we find 

• HLV: 16.4 deg^ 

• HLVA: 4.22 deg^ (-74%) 

• HLVJA: 2.78 deg^ (-83%) 

So that if LCGT is the first IFO to be added to HLV, 
there is still another 20% to be gained at the first order, 
and 65% using first and second order, by adding AIGO 
(with Indigo the differences are much smaller, to the 
close position of LCGT and Indigo on the Earth). 

Note also that the total gain going from 3 to 5 IFOs is 
not too different using the first order errors only (—71% 
with HLVJI and -78% with HLVJA), or the first and 
second (-84% with HLVJI and -83% with HLVJA). 
On fig. [9] and [lO] we plot the sky errors together with 
the errors on the arrival time. 
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Figure 9. (Color online) Skymap for a BNS signal detected 
in HLVJA. The error on the arrival time is plotted in the 
background and its value is shown in the colorbar (in ms) 
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BNS 


BHNS 


BBH 






a[l] 


a[2] 


a[l] 


a[2] 


o\\\ 


a [2] 




M [%] 


4.47e-3 


2.75e-3 


2.75e-2 


1.79e-2 


0.21 


0.15 




f] [%] 


0.21 


0.13 


0.37 


0.21 


1.56 


0.97 


:> 


ta [ms] 


0.14 


0.11 


0.22 


0.12 


0.41 


0.26 




lat [mrad] 


11.6 


9.39 


15.2 


7.68 


20.0 


9.27 




long [mrad] 


37.3 


35.6 


42.5 


26.6 


56.9 


38.8 




Af7[deg2] 


1.37 


2.60* 


2.14 


2.93* 


3.77 


4.88* 




Net. SNR 


15.93 


14.95 


14.06 






BNS 


BHNS 


BBH 






a[l] 


a[2] 


a[l] 


a[2] 


a[l] 


a [2] 




M [%] 


4.44e-3 


2.74e-3 


2.74e-2 


1.79e-2 


0.21 


0.15 


<; 


r? [%] 


0.21 


0.14 


0.37 


0.21 


1.55 


0.96 




ta [ms] 


0.11 


0.10 


0.17 


0.10 


0.35 


0.23 




lat [mrad] 


8.19 


7.47 


11.2 


7.32 


15.0 


8.51 




long [mrad] 


37.1 


38.2 


44.0 


29.2 


58.6 


39.9 






1.02 


2.78* 


1.72 


3.17* 


3.09 


4.86* 




Net. SNR 


16.18 


15.18 


14.29 



Figure 10. (Color online) Skymap for a BNS signal detected 
in HLVJI. The error on the arrival time is plotted in the 
background and its value is shown in the colorbar (in ms) 



Table VH. Same as table [y] but using 5-IFO networks 



It is also important to use the results including the 
second order of the expansions. Such role, with the 5 
dimensional parameter is illustrated in the tables V, VI, 
and VII. If we include the second order for an average 
network SNR of 15, the solid angle uncertainty increases 
on average 162%, or between 79% and 248% with 3 
IFOs, between 80% and 576% with 4 IFOs and between 
29% and 173% with 5 IFOs, depending on the Network 
and mass bin considered. 



(b) The results of this section suggest a word of cau- 
tion in reducing the parameter space and relying only 
on the IFIM. 

If someone would use the two parameter approxima- 
tion and only the Fisher information matrix, the square 
degree uncertainty would be underestimated (on aver- 
age for the examples considered in this paper) by a fac- 
tor 6.7 . 

The results of the 2 dimensional parameter space pre- 
sented in tables |ll III and 



quantitative discussion of t 



IV) are less relevant in the 



le errors because the angu- 



lar resolution is 4.5 times smaller than the values in the 
tables for the 5 dimensional parameter space presented 
here, (notice that the increase in the uncertainties with 
the size of the parameter space was already discussed 
before for a single IFO parameter estimation scenario 
based on the IFIM, e.g. in [T^). 



V. CONCLUSIONS 

In this paper we have described the most accurate 
analytical frequentist assessment to date of the uncer- 
tainties in the estimation of physical parameters from 
gravitational waves generated by non spinning binary 
systems detected by Earth-based networks of laser in- 
terferometers. The technique we adopted to quantify 
the errors is an asymptotic expansion of the uncertain- 
ties in inverse powers of the SNR ratio, where the first 
order is the inverse Fisher information matrix, and it 
provides results which are better suited to describes low 
SNR regimes. 

We have quantified how the accuracy for the intrin- 
sic parameters depends on the network SNR, while the 
measure of the direction of arrival also depends on the 



20 



network geometry. We have compared results for 6 dif- 
ferent existing and possible global networks and quanti- 
fied the relative advantages of different proposed sites. 
In particular, the fraction of the sky where the one 
sigma angular resolution is below 2 square degrees is 
shown to increase about 2.3 times when transitioning 
from the HLV network to a possible 5 sites one (keep- 
ing the network SNR fixed). Between the proposed four 
IFOs networks, HLVA guarantees smaller errors on the 
reconstruction of the position of the source. On the 
other hand, the two five IFOs network we considered, 
HLVJA and HLVJI, seem to be have comparable per- 
formances (see VIII). 
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Appendix A: Symmetric and trace free tensors and 
Gel'fand Functions 



In this appendix we provide some useful expressions 
and derivations that are used in the text, advising the 
reader to refer to [7|, |3T] and [32] for a complete anal- 
ysis of the topic. Useful information can also be found 
in dH. 

Let us consider a Cartesian frame, call it 'WF', with 
coordinates {x, y, z). We have seen before that the wave 
tensor can be written: 



r 



while m* (S) m* is just the complex conjugate. Both 
the matrices are symmetric and trace free (STF), and is 
then possible to expand them in a base of STF tensor. 
Given a STF tensor, Q, whose Cartesian coordinates, 
Qij, are given in the WF, one can expand it in the STF 
base as follows: 



Qij Qmy^ij 



(A6) 



Wij = h^^[mimj] + hx'^[mimj] 



(Al) 



where m is complex unit vector, whose components 
in WF are m = ^(l,i,0). 
Obviously: 



^[miTTij] = -{minij + m*m*) 



^ImiTrij] = —{rriimj — m*rn*) (A3) 



(A2) 

(A3) 
(A4) 



The explicit expression of m (8) to in the WF is: 
, / 1 i 



rriimj 



i -1 
\0 



(A5) 



where we use the Einstein convention on the repeated 
indexes, with m = ±2, ±1,0. The quantities Q„i are 
called the spherical components of Qij. 

The base tensors yfj" can be put in relation with the 
usual spherical harmonics Y^^'^{(j>,9) with I = 2 using 
the following formula: 



(A7) 



where are the components of the unit radial vec- 
tor n — (sin6' cos (/), sin6' sin0, cos 6*). Eq. A7 can be 



inverted to find the explicit expressions of the base ten- 
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1 


BNS 
2 


3 


1 


BHNS 
2 


3 


1 


BBH 
2 


3 


HLVJI 


CLes[l] [%] 


49.48 


79.24 


92.39 


18.34 


60.55 


80.27 


0.0 


25.60 


48.10 


CL95[1] [%] 


0.0 


4.15 


32.53 


0.0 


0.0 


6.57 


0.0 


0.0 


0.0 


CL68[l+2] [%] 


21.11 


61.59 


76.12 


5.88 


48.10 


69.90 


0.0 


11.07 


42.56 


CLg5[l + 2] [%] 


0.0 


0.0 


3.46 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


HLVJA 


CL68[1] [%] 


62.63 


96.19 


98.61 


13.84 


69.55 


93.77 


4.84 


21.45 


59.86 


CL95[1] [%] 


0.0 


5.89 


31.83 


0.0 


1.38 


8.30 


0.0 


0.0 


0.0 


CL68[l+2] [%] 


19.72 


64.71 


80.28 


7.61 


52.60 


67.47 


0.0 


10.38 


43.94 


CL95[l+2] [%] 


0.0 


0.0 


3.46 


0.0 


0.0 


3.80 


0.0 


0.0 


0.0 


HLVJ 


CL68[1] [%] 


40.83 


66.44 


80.62 


11.76 


41.18 


57.09 


1.73 


11.76 


35.64 


CL95[1] [%] 


0.0 


5.88 


28.72 


0.0 


0.00 


6.92 


0.0 


0.0 


0.0 


CL68[l+2] [%] 


15.22 


46.02 


57.43 


6.92 


24.57 


40.83 


0.0 


6.92 


22.15 


CL95[l+2] [%] 


0.0 


0.0 


3.46 


0.0 


0.0 


2.08 


0.0 


0.0 


0.0 


HLVA 


CL68[1] [%] 


45.32 


87.54 


96.54 


11.76 


59.52 


81.66 


4.15 


11.76 


43.25 


CL95[1] [%] 


0.0 


5.54 


27.68 


0.0 


1.38 


7.96 


0.0 


0.0 


0.0 


CL68[l+2] [%] 


13.49 


56.05 


72.32 


6.92 


33.22 


55.36 


0.0 


7.61 


23.86 


CL95[l+2] [%] 


0.0 


0.0 


4.15 


0.0 


0.0 


4.15 


0.0 


0.0 


0.0 


HLVI 


CL68[1] [%] 


22.84 


63.32 


85.81 


5.88 


34.60 


60.55 


0.0 


8.30 


24.91 


CL95[1] [%] 


0.0 


0.0 


5.54 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


CL68[l+2] [%] 


0.69 


36.33 


57.44 


0.0 


17.65 


40.83 


0.0 


0.69 


9.00 


CLg5[l + 2] [%] 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


HLV 


CL68[1] [%] 


13.15 


42.56 


54.32 


1.38 


11.76 


27.68 


0.0 


4.15 


11.76 


CL95[1] [%] 


0.0 


0.0 


6.22 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


CL68[l+2] [%] 


0.0 


25.26 


41.18 


0.0 


9.69 


20.07 


0.0 


0.0 


4.15 


CL95[l+2] [%] 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



Table VIII. Fraction of sources whose confidence levels of 
68% (1 standard deviation) or 95% (2 standard deviations) 
are smaller than 1, 2 or 3 square degrees, considering only 
the first order errors, or the first and second order. 




(A8) 



It is easy to verify that the 3^s also satisfy the follow- 
ing closure relation: 



Tr[3^"y'*] 



15 . 



(A9) 



and that they satisfy: 



yn* ^ („)";y- 



(AlO) 



which we used already in eq. ( A8 1 to save space 
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Eq. |A6| can be inverted to find the spherical compo- 
nents of the tensor Q: 



(All) 



Formally we can expand m ® m and m* ® m* in this 
base: 



where the only non-zero components of Mm and P„ 
are: 



(to ® m)ij = M^yi 



(A12) 
(A13) 



Combining eqs. | A1|A2|A3|AT2| and [JTS] we get: 



(A15) 



so that one recovers the expected result from |A14l 



Win = y, 



ij — ■J' I J 



-2 / "+ 








1 I ' ^ M 1 M 

-(/i+-z/ix) i -1 +-{h+ + ih-^) I -1 







(A16) 



Let us now introduce a second frame, called EF. We 
are interested in the relation between the components 
of a STF tensor in the frame WF and in EF. 

We assume that the WF is obtained from the EF with 
a rotation, parametrized with Euler angles (a,/3,7). 
The functional dependence on a set of Euler angle is 
indicated either by writing explicitly the Euler angles 
or with an arrow pointing from the fix frame to the 
rotated one. For example: 



and 



i^(a,/3,7) 



F{EF WF) 



(A17) 



(A18) 



both indicate a function F depending on the Euler 
angles that rotate EF ro WF. 

Eq. Al can still be used, providing that the compo- 
nents of m in the EF are used. These are obviously 
related with the components in the WF: 



EF^ 



i?-i(a,/3,7) 



(A19) 



where R ""^ (a, /3, 7) is the inverse rotation 3x3 matrix. 
Moreover, we can expand (to to) and (to* (8) to*), 
which are still STF: 



(to ® 'rn)ij 
(to* (8) m*)ij 



.r2™(EF^WF)X7 (A20) 
^r„2™(EF^ WF)3;,rj' (A21) 



where we have used Tab to denote the Gel 'fund func- 
tions of rank 2, which are nothing else than the spherical 
components of (to €5 to) and (to* (g) to* ) in this rotated 
frame. We have explicitly written down the dependence 
of T±2n on the rotation matrix (i.e. on the angles that 
parametrize the rotation). 

The previous relations can be inverted, to find the 
explicit expressions of the Gel'fand functions: 



23 



eral definition, of tlie ranli i functions: 



T2,„(EF^ WF) 
T_2™(EF^ WF) 



/Stt , 



-\ im * * * 



(A22) 
(A23) 



Some of the symmetries own by the rank-2 Gel'fand 
functions are easier to prove starting from the most gen- 



Ti„(a, /3, 7) ^ e-""^Pi„(cos(/3))e— (A24) 

from which the ones we used are obtained putting 
I = 2 (we shaU drop out the apex 2, as we only use 
the r ank-2 functions so no confusion is possible). In eq. 
lA24l we have introduced the functioij3 



plM = A{l-^J)-'^^l + ^i)- 



jl—n 



-[(i-A.y-'-a+A*) 



\l+m-\ 



(A25) 



where the numerical coefficient A has the value 



A; 



(■_2y-m-n-m (I ^ m)\ {I + n)\ 
2i{l-m)\ V {I + m)\ {I - n)\ 



(A26) 



Is is easy to check that P„„j has the following sym- 
metries: 



P = P = P 

~n—m nm mn 



1 \ m+n pi 



that can be used to show that: 



-n = (-)- 



-7n~nrri * 



rj-i ^ 



(A27) 



(A28) 



Finally it can be shown (^) that there is a very sim- 
ple relation between the Gel'fand function associated 
with a rotation R and that associated with the inverse 
rotation: 



(A29) 



Using eq. |A20| and |A21| on eq. |A1[ the wave tensor 
in the EF can be written: 



where it is understood that the T±2n = T±2ri(E'F — )• 

WF) =T±2n(a,/?,7)- 

One of the advantages of having a STF tensor ex- 
panded in the STF base, other than the geometrical 
insight, is that one can take advantages of all the pro- 
prieties of the Gel'fand function. In the main text, for 
example, we need to write the detector tensor of the I- 
th IFO (which is also a STF tensor) in the wave frame. 
We can apply to that tensor the same manipulations de- 
scribed in this appendix, and obtain an expression con- 
taining: T±2n(WF — ^ IDT) where we need the angles 
that rotate the Wave frame in the I-th detector frame. 
If many detectors are present in the network, however, 
it would be better to calculat e th e position of source in 
a common frame (using eq. 2.4 the Euler angles can 

We can 
^ IDT 



be obtained), which can be the Earth frame, 
achieve that result, writing the rotation WF 
in terms of the two successive rotations WF — > EF and 
EF — ^ IDF. The behavior of the Gel'fand functions over 
two successive rotations in very simple. Let Ri and i?2 
be two rotations, then the Gel'fand function associated 
with R = R2R1 can be obtained from those associated 
with i?i and R2 as follows: 



(A30) 



^ Note that 1 - n=0 does not imply that the term within deriva- 
tive does not give contribution, but only that the derivative is 

not to be performed: J^cf ^ = ^(p) 



^n^„(i?2i?l) = T,ns{R2)Tsn{Rl) (A31) 

s = -2 

which allows us to write: 

T±2ri(WF -> IDT) = r±2„(EF ^ IDF, WF ^ EF) = 
2 

= ^±2^(EF ^ IDF) r,„(WF ^ EF) (A32) 

s = -2 
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Appendix B: Proof that the generalized antenna 
pattern go to the single ones 

In this appendix we show that the generahzed an- 
tenna patterns, eq. 2.34| and eq. |2.34| do become equal 
to eq. |2.24| and eq. |2.24 if a single IFO is present in the 
network. 

If our network is composed of the I-th detector only, 
we can identify the fide frame with the detector frame, 
and then it follows from the definition 12.21 that : 



(Bl) 



because now FF=IDF. 



For the same reason (eq. 2.171: r/ = 0. 

Let us start from X2^ eq. 2.33 in which we now set 

= = 7-f = 0: 



f 1 + cos ( 



that leads to the final expression for X2' 
'(1 + cos2 0) 



cos 2-0 cos 20 — cos Q sin 2-;/' sin 20 



+ * I ^ — ^ — ^ sin 20 cos 20 + cos Q cos 2-0 sin 20 



The antenna patterns can be obtained from X2 as 
shown in the main text: 



X2 = nMO.^) (r2,(o,o,o) + r_2. (0,0,0)) (B2) 



Where 0, 6*, -0 still indicate the Euler angles that ro- 
tate the FF (i.e. the IDF) to the WF, eq. [?T] Starting 
from the general definition of the Gel'fand functions, 
eq. A24 is a matter of a few moments proving that: 



r±2.(0,0,0) = 5: 



±2 



where, as done before, we dropped out the apex, be- 
cause Z = 2 in this whole article. 



Eq. B2 then becomes: 



2-2- 



(B3) 



Using again the definition of the Gel'fand function is 
easy to find the explicit expressions of the two terms 
above: 



-Fx 



Mx2] 



(B4) 
(B5) 



from which eqs. 2.24 and 2.25 follow immediately. 



Appendix C: Second order non-diagonal elements 
of the covariance matrix and second order bias 



In this appendix we give the explicit expression of 
the second order bias and second order non-diagonal 
elements of the covariance matrix, for a network of N 
IFOs. The V to be used are the network ones, given in 
the main text, eqs. 



3.38 • • • 3.49 For the nondiagonal 



terms of the covariance matrix we find: 



'^nmpq 



^ 1 1 ' 



1^1 



'^npni'^qzt ^^^npq^mzt ~f~ ^qz,n^mpt ~l~ '^qz,m'^npt 



-\- '^Vqp.z'^nmt H" '^'^nt.z'^nipq H" ^'^mt,z'^nqp ~l~ ^^'^mt^n'^pqz H" ^'^nt^m'^pqz H" 



~l~ ^pt,7n^nq,z ~l~ '^mq,z'^pt,n ^nq,m^pt,z ~^ ^'mq,n^pt,z ~l~ ^nq,p'^mz,t 



(Cl) 
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It is easy to show, even though somehow lenghty, that ements, eq. |3.17[ follows from this one setting i — j- 
the expression we gave in the text for the diagonal el- The second order bias has the explicit expression: 



b ^-r.[2] = [Vabcde + 4^(Sac | Sbde)'' + 8^(sde | Sabc)' + 4Ua6ce,d] 

1=1 1=1 



{2VafedVgb,c + '2-VbedfVac,g + 4:VabedVgf^c) + [VafedVgcb + 

N N N 



4 

"^VabedVgcf + "^VdbegVacf) + {2Vaed^{Sgb I S/c)^ + 4Wac/ ^{^dg \ SebY + '^Vbed^{Sac \ SgfY 



1=1 1=1 1=1 

N 

+ 2w/c6 ^(Sag | SedY) + {4:Vafe,gVdb,c + '^Vafe,cVdb,g + ^Vdbe,gVaf,c) + {2Vabe,gVcdf 
1=1 

N N 

+ ■iVdbe,gVacf + ^VabejVcdg + ^Vdge.bVacf) + (4^(sQg | Sfc)^ Ved.b + 4^(Sed | Sfc)^ Vag,b 

1=1 1=1 

N -| 

+ 4:J2{Sag\SedY ^Jcm) 



1=1 

pmap??cpdep/gpfi 



VadfiVebcVgU + ^VetcVgbi + ^VgbeVtci + SVgbtVeci + ebcV gt 



8 

+ 4:VetcVgbA + 2VgtiVeb,c + 4VgtcVeb,i + SVgbtVcea + SVgbtVci,e + SVgbeVctA + S,VcteVgb,i 
+ 'iVctiVgb,e + -iVgt.iVeb^c + 'iVeb.iVgt,c + 8Vgt,bVic.e + S>Vgt.eVic,b + 4:VbetV g ,c,i) 

+ Vdcii8VbgtVaeJ + ^VbgfVae,t + SVae,tVbgJ + 8Vae,fVbg,t + 8v a f ,bV g e ,t)] (C2) 



Appendix D: Single IFO upsilon 

In this appendix we give the explicit form of the single 
IFO terms that appear in the second order covariance 
and bias. In this appendix it will be assumed that all 
the quantities refer to one of the IFOs in the network, 
and the apex (/) to unburden the notation is dropped 
out. 

Let us recall that we have defined: 



Os ds 
(a---i\j---p) = (sa-i\sj...p) = ( QQa . . . g0i\ ggj . . . QQp) 



Va,b = 


-Vab = 


Vab,c= 


{ab 1 c) 


Vabc.d = 


{abc 1 d) 


Vabc = 


-{ab\c] 


Vab,cd — 


{ab 1 cd) 


Vabcd = 


— {ab 1 a 



and that is the Fisher matrix. 



— {abc I d) — {abd \ c) — {acd \ b) — {bed \ a) 
Using this notation, the single IFO v's can be written '^"■b ,c ,d = —{a\b) {c \ d) — —Tab Fed 

my- 
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^abc , de 


{abc 1 de) — 


^ de^ abc 






^abcd , e — 


{abed \ e) 








'^abc , d , e — 


^ debate 








'^ah . cd , e — 


^ab'^cd , e 


^cd'^ab , e 






^abcde — 


— {abed \ e) 


— {abee \ d) 


— {abde \ e) 


— {aede \ b) 




— {bede \ a) 


— {abe 1 de) 


~ {abd 1 ce) 


— {aed 1 be) 




— {bed 1 ae) 


— {abe 1 ed) 


~ (ace 1 M) 


— {bee 1 ad) 




— {ade 1 be) 


— {bee 1 ae) 


— {ede 1 &a) 
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